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PREFACE 


This  report  documents  a  study  performed  by  Systems  Research  Laboratories, 
Inc.  (SRL)t  Dayton,  Ohio,  for  the  Air  Force  Aerospace  Medical  Research 
Laboratory  (AFAMRL) ,  Human  Engineering  Division,  Optical  Countermeasure 
program.  This  work  was  performed  under  Contract  F33615-79-C-0500.  The 
Contract  Monitor  was  Mr.  Donald  McKechnie  and  the  Program  Manager  was 
Maj.  Allan  M.  Dickson.  The  SRL  Project  Manager  was  Mr.  Kaile  Bishop. 

The  author  extends  his  appreciation  to  Dr.  Carroll  Day,  Messrs.  Walter 
Summers  and  Maris  Vikmanis  of  the  Human  Engineering  Division  of  the  Air 
Force  Aerospace  Medical  Research  Laboratory,  Wright-Patterson  AFB,  for 
their  valuable  comments  and  discussions. 
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Section  I 
INTRODUCTION 


The  modeling  of  human  performance  in  an  antiaircraft  artillery  (AAA)  system 
has  been  extensively  studied  by  many  investigators  in  the  past  decade  [e.g., 
Kleinman  and  Perkins  (1974),  Phatek  et  al.  (1976),  Kou  et  al.  (1978)].  Most 
of  these  works  dealt  with  the  modeling  of  human  response  in  a  simple  tracking 
task.  In  the  event  of  interrupted  observations,  the  operator's  tracking 
performance  degrades  significantly  during  the  interruption  period  and  poses 
an  appealing  modeling  problem.  The  author  tackled  this  problem  by  degrading 
•several  observer  and  controller  gains  in  the  model  and  proved  to  be  rather 
successful  (Yu  et  al.,  1980).  Efforts  were  then  directed  to  study  human 
response  in  a  manual  tracking  and  firing  task.  In  this  task,  the  operator 
(gunner)  directly  controls  the  gun  turret  and  fires  tracer  rounds  continu¬ 
ously  toward  the  target.  The  gunner  perceives  the  tracer  ending  position 
and  continuously  adjusts  weapon  pointing  in  azimuth  and  elevation  to 
minimize  the  tracer-to-target  error.  In  this  system  mode,  the  gunner  has  to 
play  both  the  role  of  a  tracker  and  a  lead  angle  computer.  The  conventional 
tracking  task  is  greatly  complicated  by  the  inclusion  of  lead  angle  estima¬ 
tion.  Wei  (1981)  developed  an  observer  gunner  model  which  treated  the 
tracer  information  as  delayed  measurements.  The  intent  of  this  paper  is  to 
extend  the  author's  previous  work  to  consider  an  even  more  general  tracking 
and  firing  scenario,  i.e.  to  consider  a  tracking  and  firing  task  subject  to 
external  measurement  interruptions. 

The  interruptions  occur,  in  the  real  world,  through  various  electronic/ 
optical  countermeasures,  weather,  or  terrain  conditions.  In  this  study, 
extensive  manned-simulation  experiments  were  conducted  at  the  Air  Force 
Aerospace  Medical  Research  Laboratory  of  Wright-Patterson  AFB,  Ohio.  A 
typical  helicopter  operational  trajectory  was  used  in  the  experiment.  The 
trajectory  consists  of  three  phases.  During  the  first  phase,  the  target  is 
standing  still  at  certain  altitude  and  is  half  masked  by  some  terrain  con¬ 
figuration.  At  the  onset  of  the  second  phase,  the  target  pops  up  for  a  full 
unmask  flight  and  moves  horizontally.  Blanking  of  target  is  administered  in 
this  phase  only.  Blanking  durations  range  from  1.5  sec,  3.0  sec,  6.0  sec, 
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and  full  blanking.  Certain  repetition  of  blanking  duration  is  also 
included. 

The  structure  of  the  gunner  model  in  a  tracer-directed  fire  system  in  Wei 
(1981)  is  adopted  and  generalized  here.  Nonlinear  ballistic  equation  is 
used  to  compute  the  loci  of  elevation  projectiles.  The  model  consists  of  a 
reduced-order  observer,  a  linear  feedback  controller,  and  a  noise  remnant 
element.  The  remnant  function  which  lumps  all  of  the  random  effects  due  to 
measurement  noise  and  human  neuromotor  response  noise  is  assumed  to  be 
Gaussian  with  its  covariance  being  a  function  of  estimated  target  velocity 
and  acceleration. 

The  effect  of  observation  interruption  is  modeled  by  exponentially  degrading 
the  observer  gain,  the  controller  gains  pertained  to  observed  states,  and 
the  bias  term  in  covariance  function.  Model  parameters  and  time  constants 
are  Identified  separately  with  respect  to  no-blanking  and  blanking  empirical 
data  via  a  least-squares  minimization  algorithm.  The  computer  simulation  of 
the  designed  model  shows  that  the  model  predicted  tracking  and  tracer  errors 
are  in  good  agreement  with  empirical  data  over  various  blanking  conditions. 


Section  II 

AAA  TRACER-DIRECTED  FIRE  SYSTEM 


In  a  tracer-directed  fire  mode,  the  gunner  perceives  both  the  tracking 
error  and  the  tracer  error  on  a  two-dimensional  visual  display.  The  tracking 
error  is  the  difference  between  the  target  angle  6^  and  the  barrel 
pointing  angle  9g.  It  is  also  referred  to  as  "lag  angle"  later  in  this 
report.  The  tracer  error  e^  is  the  difference  between  the  target  angle  0^, 
and  the  projectile  ending  angle  9  .  In  the  simulation  experiment,  the  pro¬ 
jectile  flight  path  ended  at  the  range  of  the  target.  Each  tracer  round 
disappeared  at  this  point,  0^,  from  the  display.  The  detailed  description 
of  the  configuration  is  described  in  Wei  (1981).  We  will  briefly  summarize 
the  underlying  dynamic  system  in  this  report.  Figure  1  is  the  block  diagram 
of  an  AAA  tracer-directed  fire  system. 
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I  bilking  I 


Figure  1.  Block  Diagram  of  an  AAA  Tracer-Directed  Fire  System 


At  any  given  time,  the  target  trajectory  input  eT  is  fed  into  a  visual 
display  device  and  combined  with  the  barrel  pointing  angle  0g,  as  well  as 
the  projectile  ending  angle  0  ,  to  form  error  signals  e^  and  e^.  The  human 
operator  observes  these  error  signals  and  generates  a  control  output  u  via  a 
controller,  or  H-grip,  displacement.  The  control  signal  then  drives  the 
barrel  and  rate  control  plant  for  a  new  barrel  pointing  angle  0  .  Tracer 
round  is  fired  at  this  angle  and  passes  through  the  projectile  ballistics 
computation  to  obtain  the  projectile  ending  angle  0  .  The  task  of  the 
gunner  is  to  constantly  align  the  projectile  ending  angle  to  the  target 
angle,  i.e.  to  minimize  the  tracer  error  e2*  The  dynamics  for  the  elevation 
and  the  azimuth  firing  system  are  very  similar.  In  addition,  the  elevation 
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system  can  be  decoupled  from  the  azimuth  system.  However,  the  azimuth 
system  cannot  be  separated  from  the  elevation  system  due  to  a  coupling 
factor  cos  (9g)EE  the  measurement  equation. 

By  introducing  a  state  vector  x^t)  =  [x^Ct),  xi2(t),  xi3<t)]r,  "T"  means 

"the  transpose  of,"  with  xu(t)  =  6 iT < c )  -  xi2(t)  *  0iT(t)  "  6iP(t) 

and  x^3  ■  eiT^t^’  i  =  1»  2*,  the  following  system  and  measurement  equations 

which  represent  the  underlying  tracer-directed  fire  system  can  be  derived, 
see  Wei  (1981) . 

xi  “  +  ui(t)  +  E^t)  ui(t-t)  +  Ft  0 iT ( t )  +  Gt(t)  (1) 

and 

^(t)  =  C^t)  x^t)  i  =  1,2  (2) 

where 


*If  not  otherwise  specified,  the  first  subscript  index  i  represents  the 
elevation  (i  *  1)  or  azimuth  axis  (i  =  2),  while  the  second  index  repre¬ 
sents  the  i-th  element  or  row  of  a  matrix. 


c2  *  cos  e1B(t) 

e1(t)  *  -1.34  x  (1-t)  x  [1+  (0.0052t+0.000486t2)  sin  eiB(t-x)] 
e2(t)  =  -1.28  x  (1-t) 

g^(t)  ■  (0. 0052+0. 000972t)  x  t  x  cos  0^g(t-T) 
g2(t)  =  0 


;T’  ui’  ^ii*  and  yi2  <*enote  t^ie  elevation  or  azimuth  components  of  the 
target  acceleration,  the  gunner's  control  output  and  the  observed  tracking 
error  (lag  angle)  and  tracer  error,  respectively. 


If  we  introduce  a  transformation  on  the  states  x^  and  x^2  by  x^  =  c^x^, 
x^2  *  ci  xi2>  x^  =  x^  then  Equations  (l)-(2)  can  be  rewritten  as  follows. 


A|(t)  3L[  +  uf(t)  +  E^(t)  u4  (t-r) 


+  F, 


e1T(t) 


+  G' (t) 


2^(0  -  C'  x’(t) 


< 


where 


and  x|(t)  =  jxj^t),  x^2(t),  x^(t)JT.  Notice  that  the  coupling  factor 
cos  01B  is  removed  from  the  measurement  equation  and  absorbed  into  the 
system  equation.  Equations  (3) -(4)  represent  a  nonhomogeneous  linear 
time-varying  system  with  a  time-varying  delay  in  the  control. 


Section  III 

AN  AAA  GUNNER  BLANKING  MODEL 


In  Wei  (1981),  the  author  proposed  an  observer  gunner  model  for  gunner 
performance  in  a  tracer-directed  fire  system.  The  function  of  the  gunner 
can  be  decomposed  into  two  parts  to  be  modeled.  In  the  first  part,  the 
gunner  observes  continuous  signals  and  makes  an  estimate  of  system  states 
based  on  his  internal  model  of  target  motion.  In  the  second  part,  the 
gunner  utilizes  the  observed  and  estimated  states  to  form  and  exercise  a 
control  action  in  order  to  achieve  his  objective.  The  former  one  corre¬ 
sponds  to  an  estimation  process,  while  the  latter  corresponds  to  a  control 
process.  The  reduced-order  observer  is  used  in  conjunction  with  a  linear 
feedback  control  law  to  model  the  gunner's  function.  The  structure  of  the 
model  is  shown  in  Figure  2. 


Figure  2.  Block  Diagram  of  an  AAA  Gunner  Model 

This  model  structure  is  retained  for  the  blanking  case  except  the  ballistic 
equation  is  no  longer  parameterized  by  a  linear  equation  relating  0p  and  0^. 
Instead,  a  more  realistic  nonlinear  ballistic  equation  is  used  as  shown  in 
Equations  (3)-(4).  In  addition,  time-varying  gains  are  used  to  model  the 
effect  of  observation  interruption.  We  will  discuss  the  no-blanking  case 
first,  then  the  blanking  case. 


NO  BLANKING 


The  equation  representing  the  gunner's  internal  model  of  the  tracking  and 
firing  system  can  be  written  as 

i^(t)  -  A^(t)x^(t)  +  B^(t)ui(t)  +  E^(t)  ui(t-x)  +  Gj(t)  (5) 

^(t)  -  x^<t)  i  =  1,2  (6) 

Since  both  xj^  and  xj2  are  measurable,  the  only  state  that  needs  to  be 
estimated  is  x^.  The  state  reconstructor  equation  for  can  be  derived 
by  applying  the  reduced-order  observer  theory  (Luenberger,  1971) 


^  a  •  • 

Xi3(t)  »  -(k11+k12)cix13(t)  +  k11y11<t)  +  ki2y±2(t) 

-bikilciui(t)  "  biki2ciei^t^ui^t-T^  "  ki2cigi^ 

*  kilCiCi"lyil(t)  -  kilCiCi’lyi2(t)  (7) 

The  objective  of  the  gunner  Is  to  minimize  the  tracer  error  so  that  a  maximum 
probability  of  hit  could  result.  In  other  words,  the  gunner's  response  in 
the  control  process  would  be  to  stabilize  the  underlying  system,  especially 
the  tracer  error  x^2(t);  therefore,  a  linear  feedback  control  law  of  the 
following  form  is  designed  to  achieve  this  objective. 

u^(t)  -  x^  (t)  +  vA(t)  (8) 

where 


[Y 


il’  Ti2’ 


Yi3 1 


is  a  vector  of  controller  gains  to  be  identified, 


x^(t)  -  [y^Ct),  yi2(t),  xi3(t)]T 

Is  a  vector  of  measurable  states  and  estimated  state,  v^(t)  is  a  remnant 
noise  function  assumed  to  be  Gaussian  with  zero  mean  and  a  covariance 
function 


r 

" 

[vi  (t)  vi  <s)J  ■  Kl  +  ai2 

<c) 

+  ai3 

|V'> 

6(t-s) 

(9) 

for  all  t  and  s.  are  nonnegative  model  parameters  to  be  determined. 

0iT  and  0iT  are  estimated  target  angle  rate  and  acceleration,  respectively. 
Equations  (7)  and  (8)  represent  the  gunner's  response  In  the  estimation  and 
control  process  of  the  tracking  and  firing  task.  If  we  define  a  new  state 
vector 

x^,(t)  ■  jViX Cfc) ,  y12(t),  x^^(t),  Xj^(t)  - 

then  the  state  equation  of  the  closed-loop  system  is  obtained  by  combining 
Equations  (7)  and  (8)  with  Equations  (3)  and  (4)  of  the  actual  tracking  and 
firing  system. 

^(t)  -  A±(t)X±Ct)  +  D^t)  ^  (t-T)  +  1^(0  +  E^  (t)vi(t) 

+  E11(t)vi(t-x)  +  R±(t)  dO) 

where 


The  duty  cycle  is  defined  as  the  ratio  of  the  blanking  duration  to  the  cycle 
time.  The  blanking  duration  is  the  length  of  time  that  the  target  is 
blanked  so  that  the  subject  cannot  see  the  target.  The  blanking  always 
occurs  at  the  last  portion  of  a  cycle  and  may  reoccur  periodically  over  the 
entire  TOW  firing  period.  An  example  of  blanking  sequence  is  given  in 
Figure  3. 

BLANKING 

PERIODS 


4.0  8.0  12.0  16.0  20.0  T(SFC) 

Figure  3.  Sequence  of  Blanking  for  Condition  4 

The  gunner's  performance  deteriorates  considerably  under  observation 
interruption  via  blanking  the  target.  In  Yu  (1981),  the  effect  of  blanking 
on  the  gunner's  tracking  performance  was  modeled  successfully  by  degrading 
the  gunner's  estimation  gain  k(t)  and  the  controller  gain  y(t).  A  similar 
approach  is  adopted  here  to  model  the  effect  of  blanking  in  a  more  complex 
firing  task.  More  specifically,  the  observer  gain  k^  and  controller  gains 
Yjj  and  which  pertain  to  the  observed  states  xj^  and  x^  are  assumed  to 
decrease  exponentially  as  the  blanking  starts  and  to  increase  exponentially 
as  the  blanking  stops  (see  Figure  4) . 
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During  the  recovery  period:* 


ki(C)  “  ki^cl^  + 

,ki(t0)  ‘  kl(tl)] 

/  t-t.  \ 

ail(t)  “  ail(tl> 

1  -  exp  -  J 

\  /  J 

1  - 


(15) 


(16) 


The  time  constants  x  associated  with  each  gain  parameter,  an<* 
«^(t^)  are  determined  from  the  empirical  tracking  data  collected  in  the 
blanking  experiments,  as  shown  in  the  next  section. 


*In  the  simulation  program,  the  length  of  the  recovery  period  is  defined  as 
the  minimum  of  1.5  sec  and  one-third  of  blanking  period  to  avoid  covariance 
being  negative. 


Section  IV 

PARAMETER  IDENTIFICATION  AND  SIMULATION 


The  least-squares  identification  program  developed  in  Wei  (1981)  was 
modified  to  identify  the  no-blanking  parameters.  Equation  (10)  can  be  first 
decoupled  and  then  approximated,  via  the  Average  Approximation  Method,  by 
the  following  ordinary  differential  equation  (Banks  and  Burns,  1978). 


W±(t)  -  N^t)  WjU)  +  ^(t) 


xi3(t)  = 


x14(t)  -  -k±(t)c1  xiA(t)  +  e1T(t) 


(17) 

(18) 
(19) 


where 


cici  +  biciYil(t) 


biciYi2(t) 


-1 

CiCi 


CiYil(t)ei  CiYi2(t)ei 


1 

0 

1 

T 

T 

0 

1 

0 

T 

'l  0 

’xil(t) 

0  1 

Xi2(t) 

0  0 

Si  * 

x^1(t-x) 

0  0  . 

X,12(t-T) 

1 

T 
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(l  +  biYi3)  c±  xi3(t)  -  biCi7l3xiA(t)  +  bic1vi(t) 


% 


CiXi3(t) 


Ciei(t) 


Yi3xi3(t~T)  - 


Y13xi4(t-T) 


vi ( t-T ) | 


Ci8i(t> 


IDENTIFICATION  OF  MODEL  PARAMETERS 


The  equation  which  governs  the  mean  of  states  is  obtained  by  taking 
expectation  of  Equation  (17): 


W^t)  -  N±(t)  W^t)  +  Mt  i±(t) 


(20) 


where 


wA(t) 


(1+Vi3)ciXi3(t)  '  biciYi3Xi4 


$!<*>  = 


Clx13(t) 


ciei(t)Yi3 


xi3(t~T) 


Xi4(t-T)} 


ci8i(t) 


The  first  and  second  component  of  represent  the  model  prediction  of 
ensembled  mean  of  tracking  and  tracer  error,  respectively.  On  the  other 
hand,  the  covariance  matrix  l?^(t)  satisfies  the  following  equation: 


P^t)  -  N^t)  P^t)  +  P^t)  N^U)  +  ^(t)  ^(t)  L^U) 

where 

^(t)  -  E 


(21) 


21 


! 


1^(0 


bici 


0 

0 


ciei(t) 


0 

0 


Vt} 


ail(t)  +  ai2 


- 

- 

V') 

+  ai3 

6iT(t) 

- 

0 

ail(t-T)  +  ai2 

9iT(t-T) 

+  ai3 

9  ix(t-T) 


The  first  and  second  diagonal  element  P111(t)  and  pi22(t)  of  Lj/t)  represent 
the  square  of  the  model  prediction  of  standard  deviation  of  tracking  and 
tracer  error,  respectively.  Notice  that  time-varying  parameters  are  assumed 
for  kj,  Yilt  y12»  a«d  to  reflect  the  effect  of  blanking.  Since  the 
blanking  effect  to  the  estimation  of  target  velocity  xi3(t)  is  predominantly 
expressed  through  degradation  of  k^,  there  is  no  need  to  consider  a  time- 
varying  y13,  <*12»  and  ai3> 

The  steady-state  value  of  the  parameters  are  first  identified  via  a  least- 
squares  curve-fitting  identification  program.  The  reference  curves  to  be 
fitted  are  obtained  from  empirical  tracking  and  tracer  data  collected  in  the 
manned  simulation  experiments  without  observation  interruption.  These 
experiments  were  conducted  on  an  AAA  simulator  at  the  Air  Force  Aerospace 
Medical  Research  Laboratory.  Three  simulated  helicopter  trajectories 
ranged  1500  M,  2000  M,  and  2500  M  from  the  AAA  system  were  used  as  target 
trajectories.  Figure  5  shows  some  characteristics  for  the  1500  M  trajec¬ 
tory.  Let  x^(t)  and  x^2(t)  be  tbe  ensemble  means  of  tracking 

i  and  tracer  errors  and  s^(t)  and  s^2(t)  be  the  corresponding  standard 

I  deviations.  These  empirical  means  and  standard  deviations  were  obtained  by 

averaging  and  computing  the  variance  of  the  empirical  data  from  40  simula¬ 
tion  runs  with  the  same  target  trajectory  and  the  same  subject. 

t 

i 

f 

L 
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The  parameters  were  identified  by  minimizing  the  cost  function 

defined  as  follows: 


i  =  1,2 


where  tg  is  the  initial  time  when  a  selected  tracer  round  reaches  the  range 
of  the  target,  tf  is  the  time  when  the  last  tracer  round  is  fired,  is  a 
positive  weighting  factor  chosen  to  be  one  in  the  identification  runs. 

The  direct  search  algorithm  developed  in  Wei  (1981)  was  modified  to  identify 
the  steady-state  values  of  the  parameters.  The  tracking  and  tracer  data  of 
the  helicopter  trajectory,  ranged  1500  M  from  the  AAA  simulator,  without 
blanking,  were  used  to  obtain  the  following  steady-state  parameter  values 
shown  in  Table  2. 

The  time  constants  associated  with  the  parameters  were  determined  empirically 
from  the  data  of  the  1500  M  trajectory  with  blanking  condition  5  and  listed 
in  Table  3. 


TABLE  2.  STEADY-STATE  PARAMETER  VALUES 


*Insensitive  for  elevation  case,  no  degradation  is  necessary. 


Notice  that  the  time  constants  for  t.  ,  t  ,  and  t  are  the  same  for  both 

k‘  y2  ^ 

elevation  and  azimuth  gunner  model.  This  is  as  expected  because  the  gunner 
manipulates  the  H-grip  for  elevation  and  azimuth  tracking  indiscrlminantly 
with  respect  to  the  observation  interruption.  On  the  other  hand,  a^(tg) 
and  a^(t^)  for  the  azimuth  case  are  considerably  greater  than  that  for  the 
elevation  case.  This  reflects  the  steeper  Increase  of  uncertainty  to  the 
target's  position  along  the  azimuth  axis,  because  the  azimuth  component  of 
target  acceleration  is  much  higher  than  the  elevation  component. 

SIMULATION  RESULTS 


The  gunner  model  was  implemented  on  a  CDC  CYBER  175  computer  to  simulate  the 
man- in- the- loop  AAA  tracking  and  firing  task.  For  the  convenience  of  numeri¬ 
cal  computation.  Equations  (20)  and  (21)  are  discretized  into  the  following 
form: 


where 


W  *  4>n  W°  +  Hn 
— L  -i  -i  *-i 


5“  -  £  A  K r  4 


Srt-i  ■  ‘o  +  <n+1)c 


*  W  (t  ) 
-i  ^lvn+l; 


,n+l 

“i  “  -i^n+l^ 


4^  -  exp  [N.i(tn)  A] 


(23) 

(24) 


26 


r  • 

r 


H 


■f 


®xp [t^  (t)  •  a]  do  •  M. 

l  n  —i 


exp^(tn)  * 


i±  -k(Cn  > 


2"  ‘  Si“„> 


A  =0.06  seconds 


A  simulation  program  was  developed  which  uses  the  recursive  Equations  (23) 
and  (24)  to  simulate  a  closed-loop  AAA  tracking  and  firing  task.  Inputs  to 
the  simulation  program  are  the  time  history  of  range  and  acceleration  of  the 
target  aircraft,  the  initial  angular  position  and  velocity  of  the  target, 
the  number  of  blanking  intervals,  and  the  blanking  intervals  in  chronological 
order.  Outputs  of  the  simulation  program  are  model  predicted  mean  tracking 
error  and  its  standard  deviation. 

Simulation  results  are  shown  in  Figures  6  through  17  for  the  blanking 
conditions  3,  4,  5,  6,  7,  and  10.  The  solid  curves  in  these  figures  are  the 
empirical  data  which  are  obtained  by  averaging  the  results  of  40  experi¬ 
mental  runs.  The  dashed  curve  is  the  model  prediction  of  ensembled  mean 
and  standard  deviation. 


Figure  6  and  Figure  7  show  the  comparison  of  model  versus  empirical 
elevation  mean  and  standard  deviation,  azimuth  mean  and  standard  deviation 
of  both  tracking  errors  (lag) ,  and  tracer  errors  for  the  no-blanking 


case.  Figure  10  and  Figure  11  show  the  results  for  blanking  condition  5 
which  had  a  SO  percent,  3.0  second  blanking  occur  during  [11.01,  14.01] 
seconds . 

Of  particular  Interest  is  the  comparison  of  the  empirical  standard  deviation 
in  Figure  6  and  Figure  10.  The  effect  of  blanking  the  target  to  gunner's 
performance  is  clearly  demonstrated  by  the  sharp  increase  of  the  standard 
deviation  of  tracking  errors  during  the  blanking  period  [11.01,  14.01] 
seconds.  This  effect  is  very  well  modeled  by  degrading  selected  model 
parameters  as  Indicated  in  the  model  prediction  curve  in  Figure  8.  Similar 
agreements  between  the  empirical  data  and  the  model  prediction  can  be  found 
in  Figures  8,  9,  and  12  through  17. 

These  figures  show  that  the  designed  gunner  model  can  provide  consistent 
prediction  of  the  gunner's  empirical  tracking  data  as  well  as  the  tracer 
error  data  for  both  no-blanking  and  blanking  cases.  These  figures  also 
demonstrate  that,  for  a  given  AAA  weapon  system,  the  same  set  of  parameter 
values  and  associated  time  constants  can  be  used  to  predict  the  human 
tracking  and  tracer  errors  for  all  simulated  blanking  conditions. 

However,  due  to  the  fact  that  the  helicopter  trajectory  has  very  low 
elevation  axis  maneuvering,  these  parameters  may  only  hold  for  similar 
types  of  low  maneuvering  helicopter  trajectories.  Reidentification  of 
these  parameters  may  be  needed  for  other  high  maneuvering  trajectories. 

The  computer  execution  time  of  the  overall  simulation  for  an  18  second 
helicopter  trajectory  takes  about  5.60  cp  seconds  on  a  CDC  CYBER  175 
computer . 


Elevation — No  Blanking 


Azimuth — No  Blankipg 


Figure  8a.  Mean  and  Standard  Deviation  of  Tracking  Erri 
Elevation — 1.5  Seconds,  50  Percent  Blanking 
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Figure  9b.  Mean  and  Standard  Deviation  of  Tracer  Error- 
Azimuth — 1.5  Seconds,  50  Percent  Blanking 


Figure  10b.  Mean  and  Standard  Deviation  of  Tracer  Error- 
Elevation — 3.0  Seconds,  50  Percent  Blanking 


Figure  11a.  Mean  and  Standard  Deviation  of  Tracking  Error 
Azimuth — 3.0  Seconds,  50  Percent  Blanking 
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Azimuth — 3.0  Seconds,  50  Percent  Blanking 
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Figure  12a.  Mean  and  Standard  Deviation  of  Tracking  Error- 
Elevation-^.  0  Seconds,  50  Percent  Blanking 
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Figure  13b.  Mean  and  Standard  Deviation  of  Tracer  Error 
Azimuth — 6.0  Seconds,  50  Percent  Blanking 


Figure  15b.  Mean  and  Standard  Deviation  of  Tracer  Error 
Azimuth — 1.5  Seconds,  75  Percent  Blanking 


Figure  16a.  Mean  and  Standard  Deviati 
Elevation — 1.5  Seconds,  1 


Elevation — 1.5  Seconds,  100  Percent  Blanking 
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Section  5 
CONCLUSION 

This  report  summarizes  the  modeling  of  a  gunner's  performance  in  a  complex 
AAA  tracking  and  firing  task  under  pseudorandom  observation  interruptions. 

The  highlight  of  the  task  is  that  the  gunner  fires  tracer  rounds  without  the 
aid  of  radar  and  lead  angle  computer.  Furthermore,  the  gunner's  performance 
is  greatly  hindered  by  observation  interruptions  via  blanking  the  target  on 
an  optical  display.  A  blanking  model  is  designed  which  consists  of  a  reduced- 
order  observer,  a  linear  feedback  controller,  and  a  remnant  element. 

The  gunner's  performance  is  parameterized  by  the  controller  and  estimator 
gains,  in  addition  to  the  covariance  coefficient  of  the  remnant.  The  effect 
of  blanking  is  modeled  by  degrading  these  gains  and  coefficients  as  a  func¬ 
tion  of  blanking  duration.  An  exponential  decay  form  is  assumed  for  these 
parameters.  The  associated  time  constants  are  determined  from  empirical 
data  collected  in  the  blanking  experiment.  A  direct  search  method  is  used 
to  identify  model  parameters  systematically  while  minimizing  the  least- 
squares  error  between  the  model  output  and  the  empirical  data. 

Computer  simulation  of  the  proposed  gunner  model  show  that  the  model 
predictions  are  in  good  agreement  with  empirical  data  for  various  blanking 
patterns  using  a  typical  helicopter  trajectory.  These  results  demonstrate 
that  the  model  can  adequately  describe  the  gunner's  tracking  and  firing 
characteristics  in  an  AAA  weapon  system  subject  to  observation  blanking. 

This  gunner  (blanking)  model  has  been  incorporated  into  the  MTQ  series  of 
POOl/OBS  AAA  engagement  models  and  is  designated  as  program  P001/0BS  3/6B. 

This  composite  program  POOl/OBS  3/6B  can  be  used  in  the  evaluation  of  air¬ 
craft  survivability  and  performing  weapons  effectiveness  studies.  Documen¬ 
tation  of  POOl/OBS  3/6B  is  in  preparation  at  the  Air  Force  Aerospace  Medical 
Research  Laboratory  and  will  be  distributed  separately. 


APPENDIX  A 

LISTING  OF  PARAMETER  IDENTIFICATION 
PROGRAM  ELEVATION  CASE 


•  *  u  <iu  u  «  krtucijincii 

COMMENT.  •*0«16ELIO,ID»17&OZ95,CY»H  JOE  NT  EL  PM*  USE  Hill  TR»J* 


rl  N. 

ATTACH,  TAPil,  IML6HELSUBJJ  J,CY»1,MR«1. 


DIMENSION  ALPHA  (7t ,  PS  IT  71  ,  A  (7)  ,  0(7)  •  ALP  HAN(  71 

0tMfF(JSfJN_SUM5<AriTir)7y<7,7T;r  (7,7),XI  (7,7>,3(7,7>_,F|M,:rZ(  6,TJ 
CCNHON/ARRAY/X  (1000),S(1000),ELDD  (1  b  0  0 )  ,TX«  10  00)  ,TS(  *000)  ,  A  AM  1 00  0 


CCMMON/S/C 0 ,  CEL.NSTP,  NOIN  ,T  0 ,  IPT,  Y10,  Y2  S 


INTEGER  Q, SUEXT.FR 


L1»*H 

OC  1  !•*, 7 


READ',  M.NDIP.TO,  IPT,  Y1Q0.Y10 

“pRrNr-viirKT;  ndi«;to  , iptttu - 

FCRHAT  (1H1,“NO.OF  PTS  ■  “  ,  I<*,  ZX  .“ORDER*  “,XZ» 

"C~  JETTS  I  rmi^M  l  Z75//1X  R  t  ADE  V  ER  ",  r  ?,*»  0 1  NT  Y10=-  ",  51Z 

C.S) _ 


OEL*9.03*IPT 


K*K1-KT 

Kixrrrrrr 

C  1*1,  XT) 


L^rZMNrOUMZt  AZSD,  OU  M3, 


READ  (1,  ,3)  T  ,CUH6«  AZ,  AZO,  AZOO.EL.ELO, 0UP*,AZNN.0UHZ,AZSD,DJM3 


TXCn)«XUM4 
TS(  IS  TOURS 
COMTIMUt 
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$6 


CtSTINU- 

CO*t*  34 

TrsTwpro - 

READS  (ALPHA  (It  .  1*1,7 ) 


OLOJaAJ 


FCTMAT(-OALPHAa  “,G12.5,-l“*Gtt.5,"l“,S  12. 612.  5, “l~ 


•  3*  »  IDU)3|  1  » 

PRINT  41.AJ,  IT.SU6IT 

"FOTHtn-  -Ji  ",  G 1 2.  5  *~  'ITEPAT I  CWS*TI5«''"SLBrTT’T*  TI OXS*^,  1 5) 
00  100  <*1*7 


Mill  irVUl:!  1 0 1 


CO  12  1*1*7 


UMUiCfl 


lFIU.GE.  21  .AN0.IL.LE.4I)  GO  TO  12 

rtr*=ACPHW 

CONTINUE 

CALL  IHTGIALFHANrm - 

IF  (AJ.GT.OLCJ)  GO  TO  20 


PRINT  40, ALP  HAN 


U«I*C«I  ♦£<  M 


■  JciriLafiii 


CC  15  M*l,7 


C  NTlNUc 


i wi  —  lie  ta 


GC  TO  25 


IF  IKK). LE. 1.51  A  (KJ  *0.0 


CO  33  L*l*7 


tilling 


rrviii'Hi'vi  i 


C  K*  1. 0 


IF  (CK.  (E.Q.  0)60  TO  100 


|.'i »  11  kllKFi  ( 


00  32  N*l,7 

“51MW(I*XX  (lim**2~ 
iU'12*SVM2*XI  (2,HI**2 


X1*SQ2T  (SUH1 ) 


i.rtWiLiifi'i. 


X3*X1SX2 


F0TMAT  (*1J*  SG12.57*  ALPHA*!*, G1 2.  5,*»S  512.5.*  1**612.5, 

— -  -  “ 

3*  XI (  11*  *,G  12.5/*  XI  (117X1(2)*  »,G12.5) 


MU*C.3 
F2U(M)»0.0__ 
CCSYlNU- 
CC  117  )*  1  ,J 
oo  fir  <•  i,  r 
Xi< J,K)«0.0 


DC  120  1*1, 

su<3»sujj»  t  abs  <  a  <  in > _ 

IF  (SUNT.LE.EfSIGO  TO  1000 
00  110  *«1'7 


DC  110  J*l,7 
2<1t  J>*0(NI*  U(N.J) 


ccxtinu- 

OC  140  J*  1  ?  7 


DC  140  L* 1*  7 
CC  140  <*J(7 


XI(JiLI>XIUtU«Z«tU 

cc^riNUi 


SLi4»0 .0 
OC  ISO  J*  It  7 


SUM4*SU<4fXI  (1,  Jl»*2 

siH4asa^r  isur<4i 


LO  15*  J* If  7 

V«lf  JIMltlf  JI/SUM4 


KClNr  f 
0*<0UNT-1 

D<r  i?5“  <ii70 
OC  ISO  L*  1  f 7 


F2<af  *V|IC,H>  ♦  F  2  C  Q  t  M ) 

,  170  CCMUnuJ 

FdHII.3 


cc  no  1*1.7 


1*1  ri  'I 


200  SUi5(0l*SUN5  (Q)  *8 (K OUNT, HI  2 

- ai  *3URrjjuRT(un - 

OC  215  '<■1(7 


if  j  KICK/ ii  i ; 


IF  <<DIHT.LE.7>  00  TO  159 

~TT?Tr»  1 - 

SU9IT  *  0 


*• 

00  250  <*1(7 


.  1600  CALL  EXIT 
.  fcKI 

SlflOOUTINt  1  MG  (ALPHA  f  EJ) 

. - CT51fiB  ff/WlttT/XEHP  {TOO  OTTSCTPTirOOJTEOffriOOOlTTXC  1000  J,TS  (100  or 

.  l(1Atl(1010>(T  HET  ( 100  01 


r  •€  4  C  M  l  FJ^fA  H 1  n  l  kt  I  jn  l  [ » 71  i 1  FI  I  Jfl  >  1 1  fi  tA 1 1 


0I1EMSI0N  M  t  41  f  P  <4f  4>  (P1<4,4>,PZ<4,4)  ,  ALPHA  J  7 )  (  A  (4t  6) 

1, eaor,^riYi7ebwu6irr»T^>7r*rNT{S74fimv,cQc<4i4>fOC4i  — 

2,  XT  <  10  91)  (X4  (10  00)  (EOH(10  30)  f  EO OH <  100 01 

IdUtLlHCt  (A  ( 1|  1 ) (9  fl>)  (  TEA <i(l)f  £6 ll>)  ( |£AJNF|1,1)  ^BnT  tl>) 
U#TA  HTM./ 
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c 

c  initialization 


SCAL»CO‘*2/0£L 


OC  1  J*l*  NO 


lUfflPI 


CCItiJXO. 


FCI>«0. 


i  jtfiioni  mi 


P  11 .11*0.000  02  56279 


■  •jfcf^jnmunnumn 


X  3T 1) *  0  •  0  00 


ECHT II  *X3  (II 


1FTARG.  JT.-200.)  S1*EXPTARGI 


C  COMPUTE  AND  STORE  STATES  X3  AND  X4 


9b  FORMAT  (4612.  51 


K2*«-l 

X3TXlT»T3TXin>EOOfKKI  »DE[T~ 
IFCAL’HUll  .EQ.O.I  60  TO  4 


>  Iiifiisu  firifKHllllli1UtTP^  iri 


60  T0  3 


C  COMPUTE  ANO  STORE  ESTIMATED  TARGET  VELOCITY  AN3  ACCELERATION 


&CMIKKI*X3(KK)«R4IKKI 


M  1*N0*  •  J 


C  START  INTEGRATION  LOOP 


M  T*l*Y  21 

’ISETAfl - 

CC  110  <X«lt  NSTP 


*111  |f  iliat 


IFTKK.LE.M)  60  TO  100 


IFTISET.NC.il  60  TO  101 

"S PEATS  TTTir-XEMrTRKll**2*Tlir21=TX (111 »»T - 

SSKTSORTTPT  It  111 -SEMPTtl  l**2»T  SORT  TP  T2  *211  -TST1  ll*»  2 


Vi  1  Ut>t.  ♦  .*«b*  1AUI  KK>»*  21 'ij  N<THfcMKK-M>> 

ALT»AIP"U<2> 

- (HPiWTI - 

K1««U 

- SBNsr/nunnn - 

A2*1.MTAU(K1I-TAU(KK)I/0EL 

- - 

A  (li  II  »-C  OR 
A  <I*Tl **CTi*JXPhXV3J 
A  (2.NS1)  *A(  ldl*  All*  A  2 

- 11*7*3  Hm  17“2)*ALI»A? - - 

_ 00  2  I»l,  N02 _ _ 

A  IJlf  I)  *RA 

J'lTlt  JTTV^R* 

2  CONTINUE 

- CIuTTTS  R  TTRCIH7B7  OIUTEBTCBIinT  5) - - ~ 

CR1*C0*ALPHA  (hi  *AL1*A2 
^ fii*. t fl  *4Lf*^A  tk\ 

SCAL1«SSAL*<AL1*A2I*»2 

ccTT  2  ,“2  f  ATLPhSI  5  f»S  CATI - - 

C  IFUNOOCKK,  1001.  EQ.OI.  OR.  (KK.ECl.il)  PRINT  99,  T  »  M  (1 1  ♦  Y 1 0. S1E  AN, 

C  1  U&rWUTC)  liSSfi  '  ‘ 

9a  FORMAT  C5G12.  SI 
C 

C  COMPUTE  MEAN  TRACKING  ERROR 


OC  110  I*  I*  NClH 

«  B(TT2ir~ i«1,nCIm 

l/IIMDIIl  *EA  U,J>*4(J> 

,-mi  twran  ; 

„  110  CCNTINUE 

, - FUT  5Tr.'*TJ  Rm'F  X3TTKK  I  »CRt*  XI.JKRI - - 

F<21*XJ<KK|.CRS*(X<,(KK-M»-X3<KK«MI  1*0.0  Ql»<  l. -A  21*  (5 .2*0.92  2* 

. - C  TA  U  ( K  KIT^CITSl THETTIW -HTJ - - - 

.  DC  130  I*  1.  NCIM 

« - urur  j»i;z - 

C(II*C(II*EA1NTII,JI»F(JI 

,TTB - cmw- - 

.  N  (II  *D  III 

. - - - - 

.  130  CONTINUE 

T  ■ 

c  COMPUTE  ERROR  DUE  TO  MEAN  TRACKING  ERROR 

.15 - ; - 

SKEAN«SMEAN*  (H(lt-XEMPCKK>)**2*<H<21-TX  (KPI)**2 

.,15 - 

..  c  COMPUTE  COVARIANCE  MATRIX 

."C - 

CCS  (1,11*  (ALPNA(5I*ALPHA(6)*ABS<ECH(KKI  I ♦  AL»M  A  ( 7J* A3 S  (EOOH  ( KKI 1 1 

.. - r"»?cTL - * - - 

IF(«.GT.M|  CQC(2,2I*  <  ALP  HA  <5  I  *  ALPH  A  ( 61  *ABS  <  E  OH  <  KK-1 1 1 

- 1 — ♦  Al  PMA  <  ?»  •  Af  S  If  BBWniK=^n  I  •  SCAL 1' - - - 

CALL  MUlTCEAlNT.CQS.NOIM,  N1,P1,101 
„  (.AU  MULTTEI  ,P,N0IM.N1,R2,1BI 

.  OC  220  I  *  1 .  N  Cl  N 

. - crr?TnrnTreni - -  — 

P([,JUJ1(I,  J1*P2(I,JI 


EM 

SlIRDjTINE  HU.T1E,F,L*L1.H,NR) 

- DT^E^nurmTiTur.Tjnsi  ,ncii — 

CO  10  I*ltL 

. - im - 

.  OO  10  K  *1 »  L 

OC  5  J*ltLl,L 

.  TEHF»rE  hp  *rr  i  r*r  mi 

.5  II-II  ♦  1 

„ - KK»r<“-ir*c»i - 

„  H«K>*T'MP 

,,-ni - g  K<rrn*p - 

I F(  MR.  ED.  II  RETURN 


tjc"70  nnr 
00  20  <*!•(. 


tti»*i r. 

1  I*K 

VC  i5  7STTTT7T  ' 

Tt1F*TE  1P»6(  J)*EIIII 

IS 

n»n*L 

KK»<<-1>*L»I 

L2=L-1 

- re  jg-mviz - 

L  J*I»1 

- mr-m-xfOTc - 

K1»(I-1I»L»J 

- i  j  -  ii  *m - 

3 u  N  (<  1 1  *H(K3| 

- m - 

SLRPDUTINE  OSCRT(NOIM.A»OEL.EA,EAINTiNT) 

- DnrrsnN  a  ( n  »£Aiti «  eaintt  u  .c  teft  joi — 

C  SETS  EA.EXPI  A»DEL1  ,  E  A1  NT=  I  NT  tGRAL  EA  0  TO  DEL 

- RCTTTisTrmrri - 

NMNDIM»NDIM 

- immTT=i - 

C  CEF  ( NT)  »1, 

- oo'iT'tTrrirm - 

II*NT-I 

„ - nr  crerrinrpEX»coErrri»ii/FLOATTir - 

C  NT  MUST  6c  AT  LEAST  3 

„ - "CAL''  JTACTNDIWtEAINTrAt:oEFttTreoert Tit 

CO  60  L=»3,NT 

, - C#Lr-FOtrcATEAIIir,1«JIfT,HN-,EAl  11 - 

..  IFTl.EQ.NTlGC  TO  TO 

.  - E0~CALL  DTAGTNTJTN«tAXNTTEAt  rrtt*COEFlL11 - 

70  OC  09  IX*1«NK.  NOXH1 

• - rAnir»EArm*rm - 

00  CCNTINUE 

.  - EFT) - - - 

SUO’DJTINE  0IA6IN0XH,  A,»,CltC2l 

- OrtEtETON-#  t  n  ,Bttl - 

NCtH«NDIH*l 

- NK»ioTN*ironi - 

.  NHMDH-1 

. - 1 1*  I - 

XFtCl.El.l.Q)  00  TO  10 

. - or  s~  j*r»Ttn»  rom - 


00  6  !■  it  K 

E  TlTTaTTn 

Ti.ifrfonu 

*£TU9N 

TWS - 


,2.46,2,-0.017964,-9.020511 

1,  .017491,  .124433,.  42Ji8,  .  22~4i»6F-77VI797  5E-J,.  l7302£-~3 
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.  Wir.T  JOU,  3N75  TO  a  . — 1760295TWEI,  758-39  60 - 

.  MmP(ON)  . 

.  cjhhent;  "•»cvr.B"»zic,TnsrrB02q?rcr3rT_TOENr  az  p*?r-t/sr-nerr-TRa^*« - 

.  fin. 

.  "*TT!nriTrTA'nrrT'Tvrt-&KErsT3EX53rcrwrirtt7r*-r^ - 

.  A 1  TACH.TAPE2, 1NL6T<;ACcftSU3JI3,:  Y«l,  I0*L76fl295  ,MR«l. 

*  FR33RAS  O^T  (TAPEl,TAP£2, INPUT, OLTPUT) 

.  ~CHFNSITN  ALFHA(7JVPSr(7l  ,A(7T7017)  ,ALrWfM71 - — 

oil:  Si  I  IN  SUR5(b),E(7),i)(7,7)  ,Z  (7,7)  ,XI  l7,7),3(7,7),F(6),-  2(6,7) 
_tC(H3N/APcA)  /"XMOOTJT  iTSriOOO)7A20Dt  10 00T  ,  ELS  (  TO  OOV.R  AN  ( 100  0) —  '  ~ 
l.Txaaai) ,  rs  aoooi 

CC1H3N/S/C0,  lEL  «  NST  P»  NDIN^  T07ZP71 Y1 0,  Y2  0,  AZ30 


lCGISAl  tM.PJh 
'INTEGER  T, SUEIT.FR 
LG*0 

-NFAP-rr - 

NFl?t*NJAR-l 
~l> C  1  I«l,  NPAF 


CC  1  J*1,NPAS 

KT,  J)s<I  (li  J)*fl.D  . 

IF  (I.E3.  II  V(I,J)=XI  (ItJ)Bl.D 

- zrrrJTPT.  d 

bc,j)*).: 

t  — c;stinu:  - -  - — - - — ■ - - — -  

F  R*  1 

- REJr,  i,  <1,  NOT 7 ,T D , TP T,7 1 D D ,Yrt) - - - - 

PRINT  4?, <1, NDIH.TO ,IPT,Y10 

-52 - FaMTTUHi,-xt);or  pts~»  ir, 2x-,"ORoeR«-~rr?i - 

C  2<,“INIT  TIFE*  “,G12,5//1X,  ”RE  AO  EVERY  “,I2,”  »OINT  Yl  0*  “,512 

'  "  C.51  '  -  “  - 

M«T07a.03 

K*<1-<T  - - - -  ■  •  ~  ' 

OtL»1.3T»IPT 

- TT£T*C - 

HEAD! 1,  .3) II  ,LU«4,AZ,AZD,  OUM1,  EL,  ELD,  EL  CO,  3  JM2,  ELHN,  DUM3,  •  LSD, 

- C  1*1,  XT)  -  . . .  . .  -  . .  '  '  ■  -  •  - 

uC  >4  1 1 1  ,  K 

- RE3CC1.  .3)  T  ,CUH4,  A  Z,  AZO.OUMiyf  L»Et3» ELTH)V3'JSZVr LMN,  DUM3,  E L3 D 

I F ( : 3-  (l)  )  44,46 

Tf3 - FCRMU  (12G12.6) - - 

45  iF(>n;(I-l,IFT).NE.O)  GO  TO  44 

Ih*  ( 1-1)  7IPT  «1  - - -  -- 

IF(  H.N-:.l)  GO  TO  .g 

A  Z3  0  *  6  Z  ) 


Tpi - CONTINUE - 

A  30(IHI*0UM1 

. X  (IH)  *3  JM2  • '  - - - - - 

tl'.l  Im)*El-ELMN 

- “S  (IN)  *3  )43  - - - - — - - - - - 

«  AN  (IP)  *7.5 

- 1  FT3UH4.lt  ;2877VT -RAN  (IS)  xOW147-<930-r-Tl-W)JS4) 

.  I1«RAN<  H)/ 3EL 

, - -  IFtn.LE.H)  GO  TO  44  - - - — - 

ii”T  *IS£T ♦  ! 


IF(I3£T.N£,1 )  GO  TO  44  -  - 

AZT«  l*  A! 

.  VS - CONTINUT- - - - 

K  H  3  (  2 ,  *6  >  (T  IME.ClUHS.TEL,  0UH7  ,T  ELSO,  I O*  l.KT) 

LC  A’  1*1, K  - 

KEi:  (2.  *6)  T  IME ,  0UH6,  TEL,  OUN7.TELSD 

-  -  -  I F  (  c3 r  (?l  )  48,48  - - 

4b  F  C»MH  (  ,J  12  .  5) 
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. - KSrPiTtf-T - 

ktlD»,  ULPHA  (II.I-liNPA*) 

.  RUB*7E»S 

. _  PHNT»,JPS.  <ALPM<II.J«1.NPAR> _ 

.  A  j*  0  •  0 

„  ID  DC  a  I»lt  NPAR 

„  tvrri.i 

£111*0.3 
-  A<t»i?.TJ 

.  •  continue 

.  It*SU8lT«(f 

V  29*  AZT  AU*(  A  20*7100 1 

. - £MX'T*TG  (TCPhST*  JT - 

CL1J.AJ 

.  print  4 i,TAi"f hjt fhp tts p* i .  npar) 

..40  FCRNAT  (“OALPNA*  612 .5, “l“. Gl*.5. “I*.  G  12 .5.  "I G12.  5. "l“. 

..  £  Cl  2. 5 , 'I “,~G  1275, "I  “,G  12 . 51  ”  “  - 

PRINT  41.  AJtlTtSOBir 

_  "41 - rawrrT*  *,"6i27S7*rTr nut rD?^"“7r5rwsTniri srrrr ots«"-,t5i 

,  11  OO  100  <*  1.  NPAR 

5L1ITTSJ9I  TTT 
DC  1?  L  *  1 .  NP  AR 

- IIBh IK  TO  «  Al  Fnn CT  *r  HH  *7TKT,TT~~  - ” 

IF((L.G-:.2».AN0.a.LE.4)>  GO  70  12 

- IFIAlPRTN  TlTTO  70  TOT  AlPWA  NTC1  —  ALTO  AH  a  1 - 

12  CCNTIKU- 

dUTTN  TG  ULFfcAATTM  I 
IF  CAJ.JT.OLCJ)  GO  TO  20 

- Ctl  jl*  A  r —  - - - - 

PRINT  43.  CAL  PhANlNP )  •  MP*1 .  NPARI 

PTiwr"4i;  aj,  it.sustt 

£i<>«Dt<)  *6  (  K) 

- E(< )  *  5*  E  ( Kl 

0C  15  M*1 .  NPAR 

- JTPWiTRl  •  alp  RMIK) 

is.  continue 

- irrr«i  ;cr  :t;5t  ncrsTm - 

GC  TO  25 

-J»S  EKVi*V5*E(»a  ~~ 

IF  IMKI.LE.  1.51  AOOaO.a 
■?S  CmITTO 

CiC  31  L*1,NPAR 

- mTrcTTrmn  si  co-nrro - 

CmI.O 

n - cwtrsui - 

IF  IGK.NE.O.  01  GO  TO  100 
4uTi*"SuN?*Tra 
OO  3?  >1*1, NPAR 

- srcrfsmr**  i  u  ,m  *»? - 

SL12*$UN2*XI  <2,  HI**  2 

~T2 - cwrrNtf’ - 

X  INSIST  IS  UNI  I 

- X?»51RTISUH21 - 

X3*Xt/X» 

- FRITir~H,0'LD  jrUCPffATHPr.WP-T.  WPA~W I", TIT  XT 

33  FCIHAT  c»lj*  *,G12.5/»  ALPHA*(*,  G1 2.  5.*l ».  J 12 .5,*  l*.G  12. 5 

?*  i*7Sl?.?,  *  1*761?.  S,«l»ra2.Vi» ,  G 12. 5  ,* \*r - 

3*  XICll*  *  ,G  12,5/*  XM1I/XII2)*  *,G12.5) 

6C  TO  110  “  ' 

100  CONTINUE 


OC  TO  u 
S113=0.0 


TJc'iirmTwwnir 

OC  115  1*1*  NFAR 

TttTitmr- - 

F2(L,KJ*0.0 


TT5 — cxvrrrjr 

00  117  J*  1  *  N FAR 

- DOT  1 7~  <»  1,'NFAR - 

117  XKJ,  0  =  0.0 

- - cri»(rr«i^NTJK - 

120  SU3=SU13*(AeS(0<I»H 


~if  mrTTrcE  r&  s~rgo~Tcmnnr 

OC  1  TO  N«  l.NPAR 
"XC'T30~J«  1,1.  FAR 
Z  M,  J>  =3CN>  »  VIN,  J> 

CNTINO; - 

OC  lb 0  !*  1 »  NFAR 


lbO 


~Dt  TLO  XVIJNFAR - 

OC  lbO  <* J*  NFAR 

“xttj7  n  air:  {j,L»znmrr 

ccntinux 

~S  L1L«  0. 1 - 

OC  150  J=l»  NFAR 


.  ttd — srrnr*  stttc  FXTrrrjnr»? 1 

«  SL1b=S!HT  (SU7-) 

. - DC  ~1 3F~J*T*  N  FAR - 

156  »  CL  •  JI=<I<1,  Jl/SUMb 


159 


lbO 


T7TT 


-RCJNT.2 
b*<0'JNT-l 

nDtrrrr-xvrra - 

DC  160  0*1, NFAR 

-Fni*Fn)  ♦  - xr ( koont *v (xrcr 
C04TINJ5 

“DC  170  >f«TVWFAR - 

F  21  3*  Ml  *F  (Q)*V(K,H)  ♦  F2«a.M) 

— CCTTINOi - 

Fm  =  o.a 

'CCRTIHUr - 

OC  190  1=1, NFAR 

“B  «  0  J  NT ,  r  )  ■  XTTK  C0NTTTT=F2Tarrr 
SUi5( 0) =0 • 0 

~0GT30~’1*  ITNFSR - 

SUIstOI  *SUM5  <QI»B(K0UNT,M)*»2 

- scmTDMsaRT  (suHsnn - 

to  215  1*1, NFAR 

- V  (OUNT.H)  =  £.  (KOJNT,  NT/SONFnn 

215  CC9TINU' 


150 


2ll  0 


'KCiNTAKTOTrrn - 

IF  (<3UTT.Lt.NFAR!  GO  TO  159 

~ IT*IT  *1 - 

J>UTIT  *  1 

FK*1 - - - 

CO  250  <*  1, NFAR 


Tro  *rr 
OKI  *0 . 3 
“A  KT*2.  r 


1 


I  •  ,(  lb)  i  £6  ( lb>itblNT  (16),  EM  Hi  41  ,EAI  NT  (4 ,4),  ?  (  4 >  •  CQC(  4,  4)  ,  C  (  4) 

2,  XM100J)  ,X4  (10  OOI,  EON (10  00)  .EOOHdOOOl 

“t'ffTlVlt'NCE  ■  (A  ( V,ll  ,8 1  ill  ,  IE  A(l ,  1) ,  £8  (lir,  (EtttfT  (1,1  J.EtlNT  (1)» 
oata  wt/i./ 

-~c 

.  C  INITIALIZATION  _  _ 


oC  1  1*1,  NJ _ _ _ _ _ _____  _ 

OC  1  J*l, NO 

P(t,J>*3. _ _ _ 

CKII.JItO. 

A(t,JI>].  _ 

NCIWiNSfP'/J 
NCi*vr-i 
NC»*  O? 

NAMND/E-1 

oo  n"  1*1, no 

h  (I)  *  0  . 

.  dit)*o. 

U  F  (I )  *  0  • 

-  1 1  P  (I ,  I  )  *  3  .  0  0  OlS 

I SE  T  »  0 
X3(i)*4Z00 
X4(  1)  <0. 
fc  CH  ( 1  >  *xTTii 
tC1H(l) «0. 

NCIM*ND 
N  1*  N  3 1  «» 

S1*0. 

OC  13  K<«1,NSTP 
Ki*K<*1 

LE*C03  (ELG(KK)i 
af«**gEl*  Alpha')  ii*"£9 
IF(A«.iST. -200.1  S1»EXP(ARG) 

XJ(XU*'(3(KK)*cD0(<K)*0EL 

1F(Al9H4(1) .EG.O.)  GO  TO  4 

- »4(  iriTis  i'vrscgKi  fEookkt *din - 

OC  TO  l 

4  X H(Ki)  *XL  (  KifT*ECO OTK)  «OEL 

-  3  CONTINUE 

1 FT)  <n  A  XTi  KTf- X  4(<T) 

IF(«.GE.2)ECyHCKK)*(EOM(KK)-EOH(KK-l)>/OEl 

- K=TAT(K<)70EE 

IF(«.LE.M)  GO  TO  10 

- 1  SETir  SET  n - 

IF(I3£T.NE.l >  GO  TO  10 
&U)iT  13'XO  ~  “ 

N  (»)«T23*C0 

- FUTTT*T0  mOTTO  - 

P  (»,»>*  (0.00  71507*08)  **2 

- s>*»  qs  tt<  i  it  cf-  xehp  tit)  » Fz+miz  i/cB-nrrxr)  **  z 

SSI*  (  S1RT  (P  ( 1, 1 1 l/CB-SEHP  (1J)**2MSQRT(P(2»2))  /C8-TS  ( 1))**  2 

~  .  . . 

io  ccitinue 

- TUTS - 

OC  100  <K=IST,NSTP 

Ki*«*T 
M«T»iJ(«)/0EL 
R  A*n4A/TAU  (  KKT" 

OC  2  I»l,ND2 
31*1*2“ 

A  (Jl,  I)*RA 
T(  J 1  ,J  lfSVRA 
2  CC'IT  I  SUE 

T  HSfiO*  ( ?LG(  Kl)  *ELG  ( X>0  )  /OfC 
Lt*COS(ELG(KK») 


~  1 


* 

k 
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T6«-T*lEI0»TAA(SLG(KK)  I 
SCAL«(CO»COI*‘  2/OEL 

T  «T  OTt - 

A  2»  1 .  -  ( TAU(KH -TAU(  XtO) /3 CL 
TC1R«C0»»l  PH*  (2)  *CB 
A  (It  II  ■•COR»To 


A(»,N01)«-C0«*A2 

- *T?7Nrj  ifrvyzmu - 

A(?,2)«T9 

X - 

C  COMPUTE  TRA 1SITI  CM  MATRIX  EA  AND  ITS  INTEGRAL  £AI NT 

x - 

CALL  OSCRT  (NClM,bt9EL,EB,  EBINT,  51 

- c«Twr**rPHA  ur»CB  - - 

c 

X“5T*'RT  ‘  INT  EGWTI ONTOOP - 

C 


C  1  ST»T(P(1, 1) I.SSO 

99 - mMAT(5GX2.5) - “ 

C 


CO  1?0  J*  1,  NC1H 

- orrMcm  »ca immrm - 

120  CONTINUE 

XTO - CONTINUE - - 

F  ( 1  >  a  ( C3~CR4  )•  X3(KK)  ♦CRG*  X  V  ( KKI 

- rtnAxinocrns - 

IF(«.GT.MI  F(2)aF(2)  *CR4»  <  X*.  (KK-M)  -X3(  KK-N) )  *  A? 

- 00  130  I-1PTCIH - - - 

OC  140  J» 1 1  2 

- oniaorii  *ea  iNTii,3r»rai - 

140  CONTINUE 

rmi  cm  — - - 

0(11*0. 

no — ccntj  nue - — 

c 


(T COMPUTE  XCV 1PTA NCE  KATE IX - - - 

C 

- CC3I1TD*  «CFHriST*AtT>HAr6nABSTrOH(KKn  »ALPHA  (T  >*AJS  (EODH  ( <K) ) ) 

1  •SCAl 

- CCrT2t2l«ALPHAI5l»SCAL»R2»»2 - 

IFUX.GT.M)  COC  (2*2)*  (ALPHA  (5)  4ALPHA  ( 6)  *  ABS  ( EDHf  KX”HI )  *ALP  HA  (71  • 

- 1  '1ST  ED3H  IKK-HI  J  J*SCAL*A2**? -  — ' 

CALL  MULT(EAINT,CQC,N0IM,  N1,P1*  10) 

- CALr-PUCT(EA,P,  NOIN, Nl,»2,  101 - —  ■  - 

OC  P2C  I*1,NC1N 

- CC  -220  ■)«  I ,  NCIH - 

P(t,  l)aPl(IiJ)  *P2  (I  ,  J  ) 

WT  CCITlNtT-  -  - - - — - 


)-TS  m-IST)-)- 


X'COHPUTi  ERR3»  DUETff  ST ANOARO- DEVIATION 
C 


XTMT  CONTINUE - - - 

E  DEL*  (SHEA  N*MT*SSO) 

- RETURN  . . 

EM 
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SUic  0 JT INS  NU.T  (£,F,LtLl,H,MR) 
0  HE  MS  UN  £  <  1)  ,F  (1)  iG(U)  iHIll 

octo  r*i»L 

11*1  _ 

tcm  <•*  1 »  l 

rti»=*o. 

Lir?~j*i7cm; - 

lc1P*T£lP»t(J)*F<II) 
i  i»  r  t  *  i 
KK*  ( <• 1) *  L*l 
h  {<<  I  *  TEMP 
G  CKO  *T£mP 

UTM^rn.innro3N 

GC  20  1  =  1, L 
GC  20  K*I  ,L 
T£  CP*0  • 

II*k 

00  15  J*1 , LI  ,t 

T  eM«TTcp  ,6  i  jr*rrm - 

ii*ti»l  _ _ 

KK*CK-i)»L*I 

MKO  *TSMP _ 

L  Z*L*1 

OC  31  I *1 1 L2 


L 

CO  30  J*L3,L _ 

Kl«  <I-i»*LO 
K2*»  J-1)*L»I 
H«1)«H(K2) 

in 

“i,XT»buri  n  nrairr  c  nd  r h , avb  el  ,ta,  eatnttnt  i - 

O  HEMS  I  ON  A  (l)f  EAllIf  EAINT(l)  «C  CEFC30) 

SETS  EA  =  £XMA»0£l1  ',£ArNT*tNTEGRAL  EA  BT3  OEt 
NCIM1*N1IN*1 
N  N*  N  3 1 M  *  N  0 1 H 
NH1*NT-1 
C  CS  FT  NT  I  *  1 . 

GC  10  I* 1 1  NT  PI 
II*  NT  -I 


-  10  CC£-<III=OEL*COEFCII»i)2FLOATII) 

,  C  - - NT-  MUST  6E~  AT'CEASr  3 

„  CALL  0IAG(NDIM»tAINT,A,30EF(l),  COEF  (211 

.  OC  68  t*3, NT 

CALL  KUlT  <A,EAINT,NOIM,NN,EA,  1) 

«  'lFCL.VQ.NT  I  OTTO  20 

..  bO  CAL  CHGCNUlM,EAINT,EA,1.0,COEF<U) 

..  20  GO  P  IT  It*  1,  N  K,  NOIMI 

..  £A(II)*£A(II)*1.0 

W~  ClITTFiK 

.. _ EM  _  _ 

SLI^OuTI NE'  O  IZrcrt SOIm,  *T8m7C~ZT 
0 IMFNS I  ON  A  l  1) «  S  ( 1) 

..  Net  mi*  Nit  nvi 

NMNOI  M*NOIM 

« - nRttih-i - 

.  n»i 

.  I F  (  c  r.E  1 7  IT  0  TTW  -TQ-n 

to  5  J*l, NNt NOIM 

,  GC  <•  I*  Jf  K 

- j—k  (rTiTI  lg-m - 

.  ACTI)*ACI1)»C2 

.  ^  11*  1 1 *N  )I  Ml 

. _ KETU*N _ 

10  OC  2  J*  If  NN,  NQIH 
K*J*MM1 
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APPENDIX  C 

LISTING  OF  AN  AAA  GUNNER 
MODEL  SIMULATION  PROGRAM 


M8M,T20,CN70000.  L760!95,M El, 2553960 
COMMENT. *NEM0WLSTMU6, 1 0«L 760295,  GY* 1* 

CO HMENT.  *AAA  MOPE6  3LAMKIVG  5IMJLATI0N  PROGRAM* 

AT  tACH.TAPEl,  OWL6MELSUBJ33, 13*L7  50  295,  CY*1 *NR* 1. 

FT  N, 

LGO. 

PROGRAM  S IM  J6( INPJT , OUTPJT, T APE1 ) 

COMMON/S/CO (2) ,3EL  < IH.NOIM,  V10{2),X3(2),EL«EL03,A700«MTAU,RA,A2»NO 
Al,M02,NAA.U;L,UAZ.ELTR.»2r9,I5ET ,Z1,Z2.TAU,TS(15>  ,  T  E(  15) ,T, IBL 
C  ' 

C  THE  PURPOSE  OF  THIS  PROGRAM  IS  TO  SIMULATE  AN  ELEVTN  AND  AZ INUTH  TRApK 
C  TASK  IN  THE  TRA  IE  P-3 I°ECTE3  FM£  (NODE  6)  SYSTEM 
C  SUBJECT  TO  OPTICAL  BLANKING 

C  INPUT!  THE  ELEVATION  (ED  l  AZIMUTH  ( AZ)  ANGULAR  ACCELERATION  OF 

C  TARGET.  AND  BLANKING  OURATIONSIUP  T3  15)  IN  CHRONOGICAL _ 

C  ORDER 

C  OUTPUT!  ME4N  ANO  STAND  0£V  3?  LAG  ANGLE  _ 

C  *****  ALL  ANGLES  ARE  IN  UNITS  3F  RAOIAN  ***** 

C  TAU!  DELAY  IN  SECONDS 
C  ALPHA!  PARAMETER  VECTOR 

C  ELERR!  MEAN  EL  LAG  ANGLE  (I.E.  TARGET  ANGLE-BARREL  ANGLE) 

C  AZERRl  MEAN  AZ  LAG  ANGLE 

C  ELSOI  STANOARO  DEVIATION  0*  ELEVATION  LAG  ANGLE _ 

C  AZSOI  STANOARO  DEVIATION  0?  AZ  LAG  ANGLE 

C  ELTRl  MEAN  EL  TRACER  ERROR  (TARGET  ANGLE-TRACER  ENDING  ANGLE) 

C  AZTRI  MEAN  AZ  TRACER  ERROR 

C  ELBARl  MEAN  EL  BARREL  ANGLE _ _ _ 

C  OCL  >  TIME  STEP  USED  IN  THE  INTEGRATION  ROUTINE 
C  TS( I) I  STARTING  TIME  OF  I-TH  BLANKING  DURATION 
C  TE (1)1  ENDING  TIME  IF  I-TH  BLANKING  DURATION 
C  V10I1M  INITIAL  GUESS  OF  EL  LAG  ANS.E  _ 

C  YlOZl  INITIAL  GUESS  OF  A?  LAG  ANGLE 

C  UEU  EL  CONTROL _ 

C  UAZ!  AZ  CONTROL 

C  COM.)  I  EL  RATE  CONTROL  COEFF  _  _ _ _ 

C  C0(  2)  I  AZ  RATE  CONTROL  COEFF 

C  K1 I  NO  OF  POINTS  IN  THE  ENTIRE  TRAJECTORY  _ 

C  Kl  NO  OF  POINTS  AFTER  TME  CIRST  TRACER  ROUNO  IS  FIRED 

C  ELOO!  EL  ANGULAR  ACCELERATION  3F  TARGET _ 

C  AZOOI  AZ  ANGULAR  ACCELERATION  3F  TARGET 
C  X3(l) 1  EL  ANGULAR  VELOCITY  OF  TARGET 
C  X3I 2) I  AZ  ANGULAR  VELOCITY  3F  TARGET 

C  x«»l  EST I  NATION  ERPOR  OF  ANGULAR  VELOCITY  OF  TARGET _ 

C  ELI  EL  ANGULAR  “OSITION  OF  TARGET 

C  H(l)l  MODEL  PREDICTED  LAG  ANGLE _ 

C  M( 2) I  NOOEL  PREDICTED  TRACER  ERROR 
C  P(i,lll  VARIANCE  Oe  PPE3ICTE0  LAG  ANGLE 
C  P(2,2) I  VARIANCE  0-  PRE0ICTE3  TRACER  ERROR 
C  TO!  THE  INITIAL  FIRING  TIME 
C 

REA3*.Nl,T0,IpT,  IBL _ 

PRINT  3,K1,T0,I*T 

3  FORMAT ( LM1, “NO  3F  PTS*  ", I V*  2X,*INIT  TIME*  " ,G12.5//1X,"REA0  EVERY 
C", 12,"  POINT"/) 

IF(IBL.GT.J)  READ*, (TSIO  .  TE(  K>  ,  K*l,  I9L) 

PRINT  11, IBL 

tl _ FORMAT  (IX,  I5j  IX j  "ELANKING  INTERVALS  ABE  "/) 

IFII8L.GE.1)  PRINT  i»»  ITSIO  , TE IX)  ,K*1, IBL) 

W  FORMAT  I5(1X,"(",  F9« 2,", ", ” R, 2* “)  “)  I 
KT.T0/0.03 
K»K1-KT 
T«TO 
N01M*M 
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X  PRINT  *20 /t®T 
C3(l>*1.3<* 

C0I2>*1.2« 

N01*N0IN-1 

N02*N0IM-2 

NAAsNOIM/2-1 _ 

rsr=i 

DEL=0.03»I“T 
ELSO*.005  *»3.5 
AZSO*.OQ5»»1.5 

zi-o. 

_ _ _ __ _ _ __ _ _ _ _ 

ISET*0 

UEI»0. 

UAZ-O. 

PRINT  7 

FORMAT  t/lH  .ZX.-TIMc-.SX.-EL  «E.  “  ,  9X  ,  "EL  ERR  5  X,  “E  LSD" ,  6X, 

1**EL  CTR**.  6X.  **A  7  V£L‘*«6X»‘*AZERR**.  6X."AZS0‘‘.6X.’*A2  CTR“ _ 

2.6X,~EL  TR”,6X,“A7  TR'V) 

READ (1,2) (T.OlMt, »Z, AZ3,AZOD,£L.£LO,ELOO,AZ«(N,<,ZSO,S,I*1.KT 
C  > 

00  5  1=1,  K 

READ(1,2)T,TJM1,A?,A7D,AZ33,EL.EL3.£LDD,AZMN,X,ZSD,S 

IF (EOF ( 1) ) 1, 1 _ 

IF(H00(I-1,IPT).N£.3>  33  T 3  5 
IH* ( I-i ) / I°T*1 
T*TOf<lN-l)»OEL 
TAU=7.5 

IF(0UM1.L£.287Z.  )  TAU*0UN1/ (  930.  - .  13 *OUM  1) 

NT  AU»TAU/OEL _ 

IF (1ST . EQ< 1)  DTA  0*  TAU 
IST=IST  *1  _  _ 

RA*NAA/TAU 
A2=1.-»TAU-9TA0) /OEL 
OT  AU*TAU 

X3<1)»EL0 _ 

X3(2)*AZ0 

IF((IM-1I  .SE.NTAUI  G3  TJ  3 
TAUR*TAU 

IF(0UM1.LE.  V*01>  TAUR*TAU»1NAXI(0.6,OUM1/5J30.)  _ 

YlO(l)  *-TAUR»EL3-.  001*(5.2»TAUR>.<.85*TAUR**2)*COS(  ELtO.0  5) 

Y102»-0.025»SI3S(1..AZQ) _ 

Y102*-TAUR»A73 

IF(IH.NE.l)  GO  TO  1J_  _ _ _ _ _ _ _ _ 

ELO*EL 

AZO*AZ  _ _  _ _ 

E10*Y10<1) 

_ 

GO  TO  10 

XSETaXSETH  _ 

IF(ISET.NE.l)  33  TO  10 

Zi*EL-(EL0-;l0)*.9  01»(5.2»rAUA.A#6»TAU»»2)  *SOS<  EL  0  ♦0.05) _ 

Z2=AZ-(AZ0-E20) 

_  JC  ALk_gsSEt6JLEki2Ix£WS0) _ 

ELBAR«EL-EL£RR 

IF((IM-1) .Li.NTAU)  Yt0(2)«Y102*3OS(;LfiAR) 

CALL  03SAZM  A7ERR. AZS0.EL3 AR) 

IF((N00(IH-t,I»?IYT).E7.a)  .3R.IIH.E3.D)  PRINT  6,T , LX3(  1L»_ELERR. 
l,ELSD,UEL.X3(2>,AZERR,AZS0,UAZ,iLTR,AZTR 

FORMAT (  UG12.5) _ 

CONTINUE 

FORMAT  (12G12.5)  _  _  _ 

STOP 

EMO 


SU3R0UTINE  33SE.5(EL?RR.£LS0) 

C0NN0N/S/C9  (2)  .3CL.<K.90.Y1QI2>«X30<2>  .EL.  EDO,  A  ZOO  .H.RA.AZ^ 


o  o  o 


A,N01,N02,N44,U,UA?,ELT9,  AZ  TR,  ISi T,  Z1 ,22,  TAU,  T  S  1 15  ) ,  TE  (15  1  ,T  ,  IBL 
OIN  ENSIGN  W(4) ,»(4,4>,Plt4,4),P2{4,4>, ALPHA!  7)  ,  A(  4,  4) 

1,B(16>  ,EB(16)«E3IXT!16),£A(4,4),EAI'tT(4,4),F(4)  ,CTC  (4,4)  ,  C(4) 

2«  X3(  1000)  ,  XL!  1 33  0)  ,EDH{ 10031 , ED3H( 10 00) , THET  (130  0)  ,  ALP50I1000) 
EQUIVALENCE  t  A ( 1 , 1  >  ,3(1)  )  ,  (i  A 1 1,  1)  ,  =  B(  1)  >,  (E  AINT(1 , 1)  ,EB€M  ( 1)  ) 

DATA  NT/It/ _ 

DATA  ALPHA/1.S471,  •  017491,  .  024433,  .  4 2318 1 . 22  446E-7,  .  1797SE-J, 

A  .17302E-3/ _ ; _ 

INITIALIZATION  _ _ _ _ _ 

IF(lKK-l)  ,GE.N|  S3  TO  5 _ 

if(kk.gt.ii  so  ro  5 

N1*N0**2  _ 

NOIH-ND 

ALPl-ALPi0*4LP-l»  lll_ _ _ __  _  _  _ 

ALP2-ALPHAI?) 

ALP3-ALPHAI3) _ 

ALP4»ALPH A( 4) 

ALP5-ALPHAI5)  _  _  _ 

ALP6«ALPHA(S> 

ALP7-ALPH  A(  7) _ _  _ _ _ 

ALP5Dtl)-AL®5 

ACC-0. _ 

IFLAG-0 

SCAL»C3(H»*2/3EL _ _ _ 

00  1  I-1.N0 

00  l  J-l,  NO  _ _ _ _ 

P(I, J>*0. 

CQC(  I,  J)»0. _ 

I  A ( I, J)-0. 

00  11  1=1, NO  _  _ _  _ 

Htll-O. 

0(11-0.  _ _  _ 

F(I»*0. 

II  Ptl,  I)  -0.0030 _ 

P(1,1)*0. 0010296279 

P  (2, 2)  -0.0330335677 _ _ _ _ 

X3(l  I-X30  ( II 

X4(  ll-O.  _ _ _ 

EOH(l) -X3 (1) 

EDOH ( 1) -0 . _ 

Sl-0. 

5  N(  1)  -V10(  1)  _ _ _  _ 

6  IFIISET  .EQ.l)  4(  2)  *21 

IF(IBL.LT.l)  53  TO  10  _  _ 

ALPS -ALPHA  (5) 

IS«IFLAG»1 _ 

00  12  I-IS.I3L 

IF(T «GE. TS (I).ANO.T.LT.r£(I))  G3  T3  IS 
,tt«anini(i.s,a;c/3.) 

IF(T  .GE.TE  1 1)  .4ND.T  »LT  .(TEIII+ATTII  GO  TO  16 
IF(T. GE. (TEIII»4TT>)  GO  TO  21 

_ ACC-0. _ 

IMT.LT. T$(IS)»  GO  TO  19 
GO  TO  12 

16  ALPl«ALPie*(4L<MA(l)-AL»13l  •(  1  .-£*»( -C  .4  3*  (T-TE  ( II  II) 

ALPS— 0. 0001*  (l. -EXP  (-0.43*  (T-TE(  IK  I) 

GO  TO  IS 

L  5 _ ACC-  ACC  ♦DEL _ 

IFLAG-I-1 

ALPl-ALPHA(ll*E<P(-0.0»5>-lT-TSt  I)  I) 
ALP3-ALPHA(T)*EXP(-3.52*(T-TS(I) )) 
ALPS«0.0001*C1.-EXP(-9.12MT-TS(II  II  I 
ALP 10-ALPl 
GO  IQ  46 . 


21  ALPl«AlPHACH 

ALP3«AlPNACT) 

AIP5*ALPHA<5) 

GO  TO  18 
12  CONTINUE 

18  CONTINUE _ _ 

ALP50C  KK)  «ALP5  ’ 

APG*«0EL*AL°1 

IFCARG.GT.-200.)  S1«EXPCAR3) 

tmet  cko=ei 

K1«KK*1 

K2«KK-1 _ 

C 

C  COMPUTE  TARGET  VELOCITY  ANO  ESTIMATION  ERROR 
C 

X3CK1)  aX3(<<)  ♦EOO’OEL 
IFCALP1.EQ.3.)  30  TO  4 

_  X4CK1)  *Sl*X4CKK)  ♦E0Q*C1«*51)/AL*1 _ 

60  TO  3 

4  X4CK1>*X4<«>  ♦EOO»OSL 

3  CONTINUE 

EDHC  KK) *X3  C  KK) -X  4C  KK)  _  _ 

IFCKK.GE.  21  EOON  ( KK ) * < EON(  KK) -EOHC K2 > ) / DEL 

_ X30  Cl)  «X3  CK1) _ 

C0R*C0C1)*ALP2 
C  R4»C0 ( 11  * AL°4 
IFICKK-1)  .LE.M)  GO  TO  150 

ALl*l.*0.0  01M5.2*TAU*0.48bMAy»*2)»SINCTHET  *K<r«>  » 

AC1,1)*-C0T 

AC1.2)»-G0C1)»ALP3 _ 

A  ( 2, N01) 3 A  C t,l)*ALl*A2 
AC2»ND)3A(1«2)*AL1*A2 
00  2  1*1, NO? 

J1*I *2  _  _ 

ACJ1,I)*RA 

AC Jl. J1)»~RA _ 

2  CONTINUE 

CALL  OSCRTCNO,9,OEL,EB,EBI'CT,5)  _  _ _ 

CR3*CR4*AU»»? 

SCAL13SCAL*CALI*A2)  **2  __  _ 

C 

C  COMPUTE  MEAN  TRACKING  ERROR  CI.E.  LAG  ANGLE) _ 

C 

U*AIP?*HC1>  ♦AL®3*lCC2)*ALP4»CX3(XKJ_!!«4JJCK>)___ 

00  110  I«1,N0 

00  120  J»1,'I0  _  _  _ _ _ _ 

OCZ)*QCI)  *EAC  I,  J>  *MC  J) 

120 _ CONTINUE _ 

110  CONTINUE 

F<  1  l»C  1  «“CR4)  *<3<KK)*-CR4*X4CKK> 

FC2)»X3C<K)  »CR3*  C  <4  CKK-N)  -  <3  C  KK-MI  I  »  0 . 00  1*  ( 1.- A2)  •  C  5. 2»  0.  972*f  AU 
1I«C0SCTHETCKK-D  )  _  _ 

00  130  1*1, NO 

DO  140  J*l,2 _ 

6ci)iO<I)»EAINTCI,J»*FC J) 

140  CONTINUE 
WCI)*OCI> 

OCI)«0.  _  _ _ 

130  CONTINUE 

c _ 

C  COMPUTE  COVARIANCE  MATRIX 

C  _ 

COCCI,  l)»CAt°5*AL»6*ABS(EJTCkk)l  ♦  AL»7»ABSCE0DHC  KK)  )  ) 

1  *SCAL 

COC  C2,2)«CALP*0C  K<-M)  ♦AL»5»A3S  CEO-C  C<K«M)  )  *»L  P7»  ABS  C  EOOHC  KK-N) 

1»  t*SCAll 


CALL  NULTIFAINr,  C3C, NCI, Ml,  *1,101 
CALL  MULT(*A,P,«(0,N1,P2,13( 

00  22u  1*1, SD 

00  223  J*1,V0  _  .  _ 

PII,J)*P1(I,J)*32(I,  J) 

220_  CONTINUE _ 

150  CONTINUE 

ELERR=N(i)  _ 

ELS0*SQRT (0(1,1) > 

ELTR*M(2) 

RETURN 

_ ENO _ 

SUaROUTINE  38SAT6) AZ£RR,AZ50, ELJ) 

COMMON/S/CO  (2)  ,  OEl «  <<, NO, TIO  (2>,X30(2)«EL,EL  00,  A200,H,RA  , 

1  A2  ,N01  ,ND2 «NA  1,  UEL,t),SLTR,  AZT  R,  ISET  «  21,  Z2«  TAU,  TS(  1  5)  ,TE  ( 15)  ,T,IBL 

DIMENSION  N(4)  , 3  ( 4 , 4)  *®1(4,4)  ,P2  ((*,<*)  *ALPHA(2),  A(  4 ,  4)  _ 

1,8(16)  ,£B(  15),  ESIN  T(  16)  •  £4(4, 4),  FAINT  (4,4),  F  (4)  ,COC  (4,4)  ,  0(4) 

__  2«X3(  1000)  ,X4(l)aOl  ,£0H(  1003)  .EDO HI  10  00)  ,  ALP50(in0C) _ 

EQUIVALENCE  ( »  ( 1 , 1 )  ,  3(  1) )  ,  (  EA  ( 1, 1)  ,  EB  ( 1)  ) ,  (E  AINT(  1 ,  1),EBIM(1)) 
DATA  HT/1./ 

OAT  A  ALPHA/5.5304, .11X94,.  1 27 23, 1 . 03 5 3, . 2528 6E -5,  . 19766E-3, . 757 «5E 

A- 3/  _  _  _ _ 

IF((KK-1)  .GE.H)  GO  TO  6 

_  IF(KK.GT.l)  SO  TO  S _ 

ni*no*»2 

C 

C  INITIALIZATION 
C 

00  1  1*1, NO 

_  00  1  J*1,N0 _ 

p(i,j)*'o. 

CQC  (I,  J)*0 , 

I  A (I, J) *0, 

00  11  1*1 ,N1 
N(I)*0, 

0(I)*0, 

F  ( I)  *0  , 

II  P(I,  11*0.0190 
P(1,1>*(0.'JJ59183*COS(EL3)  t  *»2 
P(2,2)*(0.3iri5Ar*COS(£L5))**2 
X3 ( 1 ) *X30  (2) 

_ X»(1)*0. _ 

EON ( 1)  *X3 ( 1) 

EOOH(l)=a. 

S1*0. 

ALP 1*ALP10  *  AL°NA ( 1 ) 

ALP2*  ALPMA(  2) 

ALP3*ALPHA(  () _ 

ALP4*ALPHA(4) 

ALP5*ALPHA(5) 

ALP6*ALPHA(6) 

ALPZ*ALPMA(7) 

ALP50(  1)  *  AL°5 

_ ACC«0. _ _ 

IFLAG*0 

C 

C  COMPUTE  ANO  STORE  STATES  X3  AND  X4 

C  _  __ 

5  M(l)*vi0(2) 

6  _ IF ( ISET.EQ.  1)  v((2)«Z2»S3S(-L6) _ 

CONTINUE 

C3*COS(ELO) 

IF(IBL.LT.l)  GO  TO  IS 
ALP5«ALPHA(3) 

IS*IFLAG*  1 
05.12  I»IS, IBL 
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ooo  cvo o  o  n  0,01  r 


IF(T.GE.TSCI)  .4'IO.T  .LT.TE(I)  )  33  T3  15 
ATT*ANINl(1.5,4;C/3.) 

IF<T.Gt.TE<I).4MD.T.LT.<rE(I)*ArT))  GO  TO  16 

IF<T.G£.  <T£<T)  *ATTH  GO  TO  21 

ACC*0. 

IFIT.LT.TSlISn  GO  TO  18 _ 

GO  TO  12 

15  ALP1*ALP10»(ALCMA(1)-ALP13IM1.-EXP(-0.43MT-TE<I>  )  )) 

ALP5*-0.001*(l.-exP(-O.43*(T-TE< 1)1)) 

GO  TO  16 

15  ACC=ACC+PEL 

_ IFLAG*I-1 _ _ 

AL>1»ALPHA«1)  *sXP«-  t.0755*«T-TS(i)  )1 
ALP2=ALPHA<2)  •  i<  P «  -  3 . 1 2*  (  T  -  T  S  (I)  I) 

ALPJ*ALP4AM)  *:<P(-0.5?*(T*TS/I))) 
ALP5*0.001*(1.-2XP«-3.12*<  r-TS«I>) )) 

ALP10=ALP1 

GO  TO  16 _ 

21  ALPl*ALpH A ( 11 

ALP2* ALPH  A  ( 21 
ALP3=ALPHA<3) 

ALP5*ALPHA<5> 

GO  TO  1? 

12 _ CONTINUE _ 

18  CONTINUE 

ALP50(KK)=A(.P5 
A  RG=-0£L*AL31 *33 

IFIASG.GT.-200.)  S1*EX°( AS J)  _  _ 

K1=KK>1 

_ K2-KK-1 _ 

X3IK11 »X3I «) ♦4TDD*0EL 

IF (  ALP1.E0.3.)  30  TO  4  _ 

XLJKllsSlMXGU'O  *A200*3jLI 
GO  TO  3 

X4(<1)»X4<«)  ♦A»PP*3Et. 

CONTINUE _ 

COMPUTE  ANO  ST02E  IS  TIMATEQ  TARGET  J/EL03ITY  AND  ACCELERATION 
E0HJK1)=X3«<1) -X4JKD 

IF (<2.G£.  1)  EtnHC<KI*(E3H(<X)-E3H<K2))/OEE 

_ X30 (  2)  *X3<<1) _ 

IFUKK-l)  .Li. Ml  GO  TO  153 

00  2  1*1 ,  N02  _  _ 

J1*I *2 

A  (  J  1«  I)  *SA  _ 

A ( Jlf Jl) *-S4 

CONTINUE  _ 

TME80* (ELG-3ELG1 /DEL 
C9*C0S (ELG) 

Tb*-TM£BO*TAN(iLG) 

SCAL*(CO(  21  *:?)  ••2/OEL 
COS*CO <  2)  *  ,L°2*3B 

_  au.1)»-co»*ts _ 

A<l,2)*-:o«2) *ALP3»C3 
A 1 2f N91 )*-C)S*4? 

AI2,N0)*A(1,2)<42 

A(2,2)«TB  . 

COMPUTE  TRANSITION  NAT  SIX  £4  4X0  ITS  INTEGRAL  E* INT _ 

CALL  OSCSTnO,l,OIL|F6,I3I'tr,5) 

c*4*co<2)*alpa»;b 

COMPUTE  MEAN  TRACKING  ESSO* 
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U*ALP2*W( II ♦  AL“(*>< ( 2)+AlP4» (X3IKK) -X4IKK)) 

00  110  I  *  1 « NO 
00  120  J^l.NO 
0(1)  *0(1)  +€A  ( I,  J )  *M (  J) 

120  CONTINUE 

11 0 _ CONTINUE _ 

FU>*(CB-CR4)*X3(KK)  ♦CR4»X4(KK) 

F(  21  *XJ(  KK)  •CB*:ru»(XM«-1)  -X3(  KK-H)  *A2 
00  130  1*1, NO 
00  140  J*  1. 2 

0(1) *0(1) ♦SAINT ( I, J)*F(J) 

140 _ CONT INUg _ 

W(I)*b(I> 

0(I)*0. 

130  CONTINUE 
C 

C  COMPUTE  COVARIANCE  MATRIX 
2 

COCTl ,  i  )*  ( AL^jTfAL  ’ 6*AB3(SH(KK> )  ♦! L  3  7 *  A  B  E7Iob H(KK)  )  ) 
1  *SCAL 

CQC  (2, 2  )*  ( ALp5 1(  K<*N)  ♦ALPB’ABS  I  SOM  (<K“M)  ) 

1  ♦  ALP7*ABS(EP0H ( < K-M) ) ) • S CAL • A2 •• 2 

CALL  NULT ( SAINT, COC, NO, Nt, *1,10) 

_ C A LL  MULT (EA .p, NO. N1.P2, 1)  I _ 

00"  220  1*1. NO 
00  220  J*1,N0 
P(I,J>*P1(I,J)432(I,J) 

220  CONTINUE 
150  CONTINUE 

_ OELG=ELG _ 

A  ZERR=H( 1 ) /CO  * 

AZSD*S')RT(»(1,1)>/CB 

AZTR*W(2)/C3 

RETURN 

ENO 

_ SUBROUTINE  NI)LT(E.F.L.Ll.H.NR) _ 

DIMENSION  Ell), r(l), 5(16), 1(1) 

00  10  1*1, L  _  _ _ 

11*1 

00  10  K*1«L  _  _  _  __ 

TEMP*0. 

00  S  J»I,Lt,L _ 

TEMP*T  EMP  ♦S(J)*r(II) 

5  __ 

KK* (K-1)»L*I 

H(KK)*TENP  _ 

10  G(KK)*TEMP 

_ IF(MR.EQ.l)  RETURN _ 

00  20  1*1 ,L 

00  20  K*I,L  _ 

TENP«0. 

II*K  _  _  _  _  _  _ 

00  15  J*I,L1,L 

_ TEMP; TEMP»G(  J)  *S  (I  I) _ 

15  II*II^L 

KK* ( K-i ) *L*t  __ 

20  H(KK)*TEMP 

(.2*1-1 

00  30  1*1 ,L2 

_ L3*J  ♦  l _ 

00  30  J*L3,L 
K1*(I-1)*L*J 
K2*( J-1)*L*I 
30  Ml K1 ) *M( K2) 

ENO 

$U}ROyTINE  OSCRT (NOIM,A,OEL,EA,EAINT,NTI 


DIMENSION  All) (-All) ,  EAU  T ( 1 ) , C3EF I  30) 

C  SETS  EA*SXsH*OEO»E4I'tr»INT£CR4L  EA  0  TO  DEL 

NOIHl  =  NOIMn 
NN=NOIi*NOIM 
NTMlsNT-1 

COEF<NT)»l. _ 

'00  10  I=1,NTM1 
II*NT-I 

10  COEFIII>*OEL»C3"F(  IU1)/*l3AT(  II 
C  ilT  MUST  3E  AT  LEAST  3 

CALL  OIAGIN3H,SAINT,A,:3£J(l)  ,:0£F(2)) 

DO  60  L-3.NT _ 

CALL  MULTIA, SAINT, NOIN, NN,  EA,1) 

IFIL.EQ.NTISO  n  70 

61  CALL  OIAGIMSIM, SAINT, EA,1.1,C0EF(L)I 
70  00  80 

EA(II)=EAIIt)*l.O 

80  CONTINUE _ 

ENO 

SU3ROUTINE  OTASI  NT  IN , A , 3 , 3 1, C2) 

DIMENSION  All), 3(11 

NDIM1=N0IM+1 

NN*NOIM»NDIM 

_ N  MlaNDIM-1 _ 

II  *  1 

IF (Cl, EQ, 1,1)  S3  TO  10 
00  5  J=1,NN,M0TN 
J+NM1 

00  4  1=J,K 

_ 4  A  1 1) »C1*9 1 II _ 

A(II)*A(II)*C2 

5  II=II*NQIM1 
RETURN 

10  00  7  J=1,NN,N0I1 
K*J*NMt 

00  6  I*  J,  K _ 

6  All) *9 (I) 

A(II)*A(II)*C? 

7  II=II*NDIM1 
RETURN 

ENO 
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